Neutrino mass bounds from cosmology 
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Cosmology is at present one of the most powerful probes of neutrino properties. The advent of precision data 
from the cosmic microwave background and large scale structure has allowed for a very strong bound on the 
neutrino mass. Here, I review the status of cosmological bounds on neutrino properties with emphasis on mass 
bounds on light neutrinos. 



Neutrinos are the among the most abundant 
particles in the universe. This means that they 
have a profound impact on many different aspects 
of cosmology, from the question of leptogenesis in 
the very early universe, over big bang nucleosyn- 
thesis, to late time structure formation. 

At late times (T < Tevv) neutrinos mainly in- 
fluence cosmology because of their energy density 
and, even later, their mass. 

The absolute value of neutrino masses are very 
difficult to measure experimentally. On the other 
hand, mass differences between neutrino mass 
eigenstates, (7711,7712,7713), can be measured in 
neutrino oscillation experiments. 

The combination of all currently available data 
suggests two important mass differences in the 
neutrino mass hierarchy. The solar mass differ- 
ence of 5m\2 — 8 X 10~^ eV^ and the atmospheric 
mass difference 5mlj, ~ 2.6 x IQ-^ eV^ |1l2l,Sl4) . 

In the simplest case where neutrino masses are 
hierarchical these results suggest that mi ^ 0, 
7772 ~ SrUsoiuLT, and 7773 ~ (5matmosphcric M- If tl^e 
hierarchy is inverted one instead finds TO3 ~ 0, 

1712 ~ I^TTlatmosphcric , and TOi ^ (^Tn-atmosphoric • 

However, it is also possible that neutrino masses 
are degenerate, mi ~ 7712 ms > (Jmatmosphoric, 
in which case oscillation experiments are not use- 
ful for determining the absolute mass scale (S]- 

Experiments which rely on kinematical effects 
of the neutrino mass offer the strongest probe of 



this overall mass scale. Tritium decay measure- 
ments have been able to put an upper limit on 
the electron neutrino mass of 2.3 eV (95% conf.) 
13 (see also contribution by G. Drexlin to the 
present volume). However, cosmology at present 
yields a much stronger limit which is also based 
on the kinematics of neutrino mass. 

Very interestingly there is also a claim of direct 
detection of neutrinoless double beta decay in the 
Heidelberg-Moscow experiment |7], correspond- 
ing to an effective neutrino mass in the 0.1 — 0.9 
eV range. If this result is confirmed then it shows 
that neutrino masses are almost degenerate and 
well within reach of cosmological detection in the 
near future. 

Neutrinos are not the only possibility for sta- 
ble eV-mass particles in the universe. There are 
numerous other candidates, such as axions and 
majorons, which might be present. As will be dis- 
cussed later, the same cosmological mass bounds 
can be applied to any generic light particle which 
was once in thermal equilibrium. 

Here I focus mainly on the issue of cosmologi- 
cal mass bounds. Much more detailed reviews of 
neutrino cosmology can for instance be found in 

IHEl. 
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only become non-relativistic around the epoch 
of recombination, significantly later than matter- 
radiation equality. While they are relativistic, 
neutrinos free-stream a distance roughly given by 
A = 27r/fc = cr, where r is conformal time. The 
total free-streaming length is therefore roughly 
given by A/ ct(T = m^). At scales smaller 
than this, the solution to the Boltzmann equation 
is exponentially damped. For scales larger than 
the free-streaming length neutrino perturbations 
are unaffected. It is therefore possible to divide 
the solution into two distinct regimes: 1) fc > 
t(T — m): Neutrino perturbations are exponen- 
tially damped 2) k < t(T — to): Neutrino pertur- 
bations follow the CDM perturbations. Calculat- 
ing the free streaming wavenumber in a flat CDM 
cosmology leads to the simple numerical relation 
(applicable only for T^q ^ to ;» Tq) 



Figure 1. The transfer function T{k,t = Iq) 
for various different neutrino masses. The solid 
(black) line is for to^ = 0, the long-dashed for 
TOjy = 0.3 eV, and the dashed for rrii, = 1 eV. 



1. Neutrinos in structure formation 

The temperature of standard model neutrinos 
IS given by = {i/liy/^T^, and the number 
density of a given flavour is therefore — 7 /8 x 
(4/11) X rij. From this the present contribution 
to the matter density from massive neutrinos is 



92.5 eV 



(1) 



where N^, is the number of neutrino species. 
Thus, even a sub-eV neutrino mass gives a signif- 
icant neutrino contribution to the energy density 
and therefore has an effect on structure forma- 
tion. 

Perturbations in the neutrino distribution can 
be followed by solving the Boltzmann equation [?] 
through the epoch of structure formation. There 
are publicly available codes, such as CMBFAST 
which can do this. 

The main difference between neutrinos and 
cold dark matter is that neutrinos are light and 



^FS 



20 Mpc /T^r 



1 + log 3.9 



(2) 



In Fig.^transfer functions for various different 
neutrino masses in a flat ACDM universe (fl^ + 
+ r^A = 1) are plotted. The parameters used 
were 0.^ = 0.04, VLcdm = 0.26 - Vl^, VLj^ = 0.7, 
h = 0.7, and n — I. 

When measuring fluctuations it is customary 
to use the power spectrum, P{k,T), defined as 



P{k,T) = |4(t)|= 



(3) 



The power spectrum can be decomposed into a 
primordial part, Po(^)i and a transfer function 

r(fc,T), 



P(fc,T) =Po(fc)T(fc,r). 



(4) 



The transfer function at a particular time is found 
by solving the Boltzmann equation for S{t). 

At scales much smaller than the free-streaming 
scale the present matter power spectrum is sup- 
pressed roughly by the factor 



AF(fc) _ Ar(fc,T To) 



P{k) 



T{k,T = To) 



(5) 
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as long as 0,^, il,„. The numerical factor 
8 is derived from a numerical solution of the 
Boltzmann equation, but the general structure of 
the equation is simple to understand. At scales 
smaller than the free-streaming scale the neutrino 
perturbations are washed out completely, leaving 
only perturbations in the non-relativistic matter 
(CDM and baryons). Therefore the relative sup- 
pression of power is proportional to the ratio of 
neutrino energy density to the overall matter den- 
sity. Clearly the above relation only applies when 
rii, r2„j, when becomes dominant the spec- 
trum suppression becomes exponential as in the 
pure hot dark matter model. This effect is shown 
for different neutrino masses in Fig. ^ 

The effect of massive neutrinos on structure 
formation only applies to the scales below the 
free-streaming length. For neutrinos with masses 
of several eV the free-streaming scale is smaller 
than the scales which can be probed using present 
CMB data and therefore the power spectrum sup- 
pression can be seen only in large scale structure 
data. On the other hand, neutrinos of sub-eV 
mass behave almost like a relativistic neutrino 
species for CMB considerations. The main ef- 
fect of a small neutrino mass on the CMB is that 
it leads to an enhanced early ISW effect. The 
reason is that the ratio of radiation to matter at 
recombination becomes larger because a sub-eV 
neutrino is still relativistic or semi-relativistic at 
recombination. With the WMAP data alone it 
is very difficult to constrain the neutrino mass, 
and to achieve a constraint which is competitive 
with current experimental bounds it is necessary 
to include LSS data from 2dF or SDSS. When 
this is done the bound becomes very strong, some- 
where in the range of 1 eV for the sum of neutrino 
masses, depending on assumptions about priors. 
This bound can be strengthened even further by 
including data from the Lyman-a forest and as- 
sumptions about bias. In this case the bound on 
the sum of neutrino masses becomes as low as 
0.4-0.6 eV. 

In Fig. 121 Ax^ is shown for an analysis which 
includes the WMA P CM B data [E], the SDSS 
galaxy survey data jl8ll4| . the Riess et al. SNI-a 
"gold" sample [TC], and the Lyman-a forest data 
from Croft et al. |17| . For the Lyman-a for- 




(eV) 



Figure 2. Ax^ as a function of the the sum of 
neutrino masses rriy. 



est data, the error bars on the last three data 
points have been increased in the same fashion 
as was done by the WMAP collaboration pj], in 
order to make them compatible with the analy- 
sis of Gnedin and Hamilton |18j . In addition to 
the neutrino mass, which I take to be distributed 
in three degenerate species, I take the minimum 
standard model with 6 parameters: fim, the mat- 
ter density, fJfc, the baryon density, iJo, the Hub- 
ble parameter, and r, the optical depth to reion- 
ization. The normalization of both CMB, LSS, 
and Ly-a spectra are taken to be free and unre- 
lated parameters. The priors used are given in 
Tabled 

In this particular analysis, the 95% C.L. upper 
bound on the sum of neutrino masses is 0.55 eV. 

In Table |21 the present upper bound on the 
neutrino mass from various analyses is quoted, 
as well as the assumptions going into the deriva- 
tion. As can be gauged from this table, a fairly 
robust bound on the sum of neutrino masses is at 
present somewhere around 0.5-1 eV, depending 
on the specific priors and data sets used. 
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Table 2 

Various recent limits on the neutrino mass from cosmology and the data sets used in deriving them. 1: 
WMAP data, 2: Other CMB data, 3: 2dF data, 4: Constraint on erg (different in 4°, 4^ and A"), 5: 
SDSS data, 6: Constraint on Hq, 7: Constraint from Lyman-a forest. 



Ref. Bound on ^ Data used 



Spergel et al. (WMAP) [H] 


0.69 eV 


l,2,3,4-,6, 7 


Hannestad 


1.01 eV 


1,2,3,6 


Allen, Smith and Bridle pH) 




l,2,3,4^6 


Tegmark et al. (SDSS) [TH 


1.8 eV 


1,5 


Barger et al. 


0.75 eV 


1,2,3,5,6 


Crotty, Lesgourgues and Pastor [22] 


1.0 (0.6) eV 


1,2,3,5 (6) 


Seljak et al. P] 


0.42 eV 


l,2,4^5,6,7 


Fogh et al. [21 


0.5 eV 


1,2,3,4^5,6,7 


The present work ESI 


0.65 eV 


1,5,6,7 



Table 1 

Priors on cosmological parameters used in the 



likelihood analysis. 


Parameter 


Prior 


Distribution 


= -\- r^x 


1 


Fixed 


h 


0.72 ±0.08 


Gaussian J^l 


nbh^ 


0.014-0.040 


Top hat 


Us 


0.6-1.4 


Top hat 


T 


0-1 


Top hat 


Q 




Free 


b 




Free 



2. General thermal relics 

While light neutrinos are the canonical light, 
thermally produced particles, there are other 
species which can have the same properties. At 
T ~ Mpi gravitons are presumably in thermal 
equilibrium. However, since inflation occurs at 
a lower temperature, there should be no thermal 
graviton background present. 

Other examples include majorons, axions, etc. 
Any such species which was once in equilibrium 
is characterized by only two quantities, its mass, 
mx, and the temperature at which it decoupled 
from thermal equilibrium, Td- At this temper- 
ature the number of degrees of freedom was g*. 
The relevant masses are <C MeV so that the parti- 
cles decouple when they are relativistic, i.e. at de- 
coupling they are characterized by a Fermi-Dirac 
or Bose-Einstein distribution of temperature Td- 



In the present-day universe, the new particles 
will be non-relativistic and contribute a matter 
fraction 

^ j{2. ^ mxgx ^ f 1 fermions 
^ 183 eV \ I bosons 

where gx is the number of the particle's inter- 
nal degrees of freedom while g^,^ is the effective 
number of thermal degrees of freedom when ordi- 
nary neutrinos freeze out with g,^ = 10.75 in the 
absence of new particles. 

In the late epochs that are important for struc- 
ture formation, the momentum distribution of the 
new particles is characterized by a thermal distri- 
bution with temperature Tx that is given by 



In Ref. 26. CMBFAST was modified to incor- 
porate such thermal relics, both fermionic and 
bosonic. 

Here I present an updated analysis for the 
fermionic case with gx ^ "2 (i.e. a single Majo- 
rana fermion), using the same data sets as de- 
scribed above. Fig. |3 shows the 68% and 95% 
likelihood contours for Qxh^ and g*- Overlayed 
are isocontours for particle masses (in eV). 

As an example, a species freezing out at the 
electroweak transition temperature has = 
106.75, and from the figure it can be seen that 
the upper bound on the mass of such a particle 
is around 5 eV. For a corresponding scalar the 
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mass bound will be slightly higher. For a species 
decoupling after the QCD phase transition where 
g.t < 20 the mass bound is roughly m<l eV. 
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Figure 3. The fuU lines show 68% and 95% confi- 
dence regions in the (r2,y/i^,g,) plane for the case 
where the thermal relic is a fermion with g — 2. 



3. Conclusion 

I have reviewed the present status of cosmo- 
logical mass bounds on neutrinos and other light, 
thermally produced particles. Already now these 
bounds are about an order of magnitude stronger 
than current laboratory limits on the neutrino 
mass, albeit more model dependent. In the com- 
ing years a wealth of new cosmological data will 
become available from such experiments as the 
Planck Surveyor CMB satellite. This is likely to 
allow for a measurement of the (sum of) neutrino 
masses in the 0.1 eV regime ^27 28 29 30 . To- 
gether with direct detection experiments like KA- 
TRIN and future neutrinoless double beta de- 
cay experiments, cosmology will provide the an- 
swer to whether neutrino masses are hierarchical. 
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